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The time evolution of tlie magnetization of a magnetic molecular crystal is obtained in an ex- 
ternal time- dependent magnetic field, with sweep rates in the kT/s range. We present the 'exact 
numerical' solution of the time dependent Schrodinger equation, and show that the steps in the hys- 
teresis curve can be described as a sequence of two-level transitions between adiabatic states. The 
multilevel nature of the problem causes the transition probabilities to deviate significantly from the 
predictions of the Landau-Zener-Stiickelberg model. These calculations allow the introduction of an 
efficient approximation method that accurately reproduces the exact results. When including phase 
relaxation by means of an appropriate master equation, we observe an interplay between coherent 
dynamics and decoherence. This decreases the size of the magnetization steps at the transitions, 
but does not modify qualitatively the physical picture obtained without relaxation. 
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I. INTRODUCTION 

There has been an increased interest in the study of 
crystals consisting of high-spin molecules such as Mni2- 
Ac and FesO. These organic molecules, also known as 
molecular nanomagnetsji contain transition metal atoms 
with strongly exchange-coupled spins, which causes the 
individual molecules to behave as a single, large spin. 
Experiments on the magnetization dynamics of these 
molecular crystals have shown the presence of a se- 
ries of steps in the hysteresis curve at sufficiently low 
temperatures . ^'^'^ This behavior is a consequence of 
quantum mechanical tunneling of spin states through the 
anisotropy energy barrier and occurs when the external 
field brings two levels at different sides of the barrier into 
resonance. This macroscopic quantum effect has been 
the subject of intensive experimental and theoretical in- 
vestigation which revealed additional remarkable prop- 
erties of the magnetic molecules. In particular, it has 
been demonstrated that the magnetic tunneling could be 
accompanied by the emission of electromagnetic radia- 
tion, and bursts of microwave pulses have been detected 
in recent experimentsi^iSii It has been proposed that the 
physical mechanism responsible for this radiation is a col- 
lective quantum effect known as superradianceiSiSiiLiLiS 
This interpretation has been questioned and it was ar- 
gued that when one includes the time scale of relax- 
ation, a maser-like effect is more likely responsible for 
the observations!"^ Furthermore, the change in magne- 
tization can be described in terms of avalanches, which 
were recently shown to propagate through the crystal in 
an analogous way to that of a flame front in a flammable 
chemical substance (deflagration) pi^ii^ It has also been 
suggestedi^ that these molecules can be used for imple- 
menting a quantum computational algorithm. 

In this work we study the dynamics of the multi- 
level system corresponding to the 21 spin states of the 
Mni2-Ac molecule {S = 10) in a time-dependent mag- 



netic field. An 'exact numerical' solution of the relevant 
time-dependent Schrodinger equation is obtained in the 
whole time interval of interest. The results show that, 
although the usual qualitative picture of consecutive two- 
level transitions holds, the transition probabilities de- 
viate significantly from the predictions of the Landau- 
Zener-Stiickelberg^7d?d? (LZS) model. This deviation 
is closely related to the multilevel nature of the problem: 
examples are given where the single relevant parameter of 
the LZS model (which is essentially the level splitting in 
appropriate units) is the same, but the transition proba- 
bilities are different due to the different time dependence 
of the adiabatic energy levels. An efficient approxima- 
tion method based on non-LZS two-level transitions is 
introduced, which is able to describe the dynamics with 
high accuracy. Furthermore, relaxation effects are also 
included. By considering realistic dephasing rates, it is 
shown that for field sweep rates in the kT/s range, nei- 
ther unitary time evolution nor relaxation dominates the 
dynamics. The interplay between these two processes re- 
sults in a decrease of the transition probability at a given 
avoided level crossing. 

The paper is organized as follows: In Sec. the rele- 
vant Hamiltonian is discussed, along with its level struc- 
ture, and the dynamical equations in the adiabatic basis 
are introduced. Results related to the solution of the time 
dependent Schrodinger equation are presented in Sec. lIIII 
and the consequences of relaxation effects are discussed 
in Sec. II VI Finally in Sec. the results are summarized 
and conclusions are presented. 



II. MAGNETIC LEVEL STRUCTURE AND 
DYNAMICAL EQUATIONS 

ExperimentaliSi2iSfli2iiS2iSi studies on crystals of 
Mni2Ac and FegO suggest that the spin Hamiltonian for 
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these systems can be written as 

Hs{t) ^ Ho{t) + Hi{t), 



(1) 



where Hq is diagonal in the eigenbasis {|rr!,)} of the z 
component of the spin operator, Sz- 



Ho{t) = -DSl - FSt - g^iBB{t)Sz. 



(2) 



Here the last term in the right-hand side describes the 
coupling to an external magnetic field applied along the 
z direction, which is parallel with the easy axis of the 
crystal. This external field is time-dependent, with sweep 
rates on the kT/s scale.^ Hi in the Hamiltonian contains 
terms^*^ that do not commute with Sz- 

Hi = CiSl + St) + E{Sl+Sl)/2 + K{S+ + S-)/2. (3) 

In the present paper we will concentrate on Mni2-Ac, 
which can be considered as a representative example of 
molecular nanomagnets. In this case the values of the 
parameters in Hq are D/ks = O-SO-fT, and F/ks = 
1.1 • \Q~^K. The coefficients in Hi, which are essen- 
tial for the determination of the transition probabilities, 
can be obtained by fitting the theoretical results to ex- 
perimental magnetization curves. ^'^ In this paper we use 
K = OmbgnsB, E/ke = -4.48 IQ-^ K, C/ke = 1.7 
10^^ K unless otherwise stated. 

Considering the total Hamiltonian as the generator 
of the time evolution, the corresponding time-dependent 
Schrodinger equation governs the dynamics. We can also 
use a density operator g to describe the system and write 



dg 
dt 



(4) 



where h = 1. Relaxation effects can then be included 
through additional terms on the right hand side of this 
equation, see Sec. llVl for more details. 

A direct calculation of the time-dependent solutions of 
Eq. Q when expanded in the eigenbasis of the spin 

operator Sz turns out to be a rather difficult problem: 
even for large field sweep rates (i.e, kT/s), the satura- 
tion of the magnetization is reached in a few milliseconds. 
During this time, roughly 10^ Bohr oscillations take place 
due to Hq, raising demanding requirements on the accu- 
racy of the numerical process. An alternative and more 
efficient way of dealing with this equation is based on the 
expansion of the time-dependent states in an adiabatic 
basis, i.e., the basis of the instantaneous eigenstates of 
Hs{t): 



Hs{t)\E^{t)) ^ E,,{t)\E^{t)) . 



(5) 



It is convenient to label these states so that they cor- 
respond to the eigenenergies in increasing order: E[j(t) < 
Ei{t) < ... < E2o{t) at all times. The energy curves 
En(t) are thus continuous. It must be stressed that the 
time dependence of these states is parametrical. In fact 



they are in a one-to-one correspondence with the exter- 
nal field, and thus the value of the time-dependent B 
completely determines the states \En). Next, the density 
operator is expanded in the time-dependent basis deter- 
mined by Eq. jSj: 



g{t) = 



M \E^{t)) {E„M ■ (6) 



Using Eq. to calculate the dynamics, and Eq. (O to 
obtain the time dependence of the adiabatic states, one 
finds that the time evolution of the matrix p takes the 
form of a von Neumann equation 



dp 
dt 



H,p 



where H is given by 



~ OHg 

Hnm{t) = i {En \ \Bn 



E„{t)~E,n{t) ' 



(7) 



(8) 



if n 7^ m, and = 0. This expression explicitly 

shows that an appreciable change in the populations 
Pnn and pmm IS cxpccted around the avoided cross- 
ing of the levels En and Em, i.e., when the denomina- 
tor in Eq. (jSJ has a local minimum. A qualitatively 
similar conclusion follows from a degenerate perturba- 
tion calculation around the avoided level 
crossings; but note that Eq. (jSJ is an exact result. 

The level scheme and the minimal energy difference be- 
tween levels at the avoided crossings are shown in Fig.^ 
As a guiding line, far from the crossings we can associate 
a label m to each energy eigenvalue En in such a way that 
the overlap |(m \ En) \ is maximal over all possible values 
of m between -10 and 10. This assignment is based on the 
fact that Hi is a relatively weak perturbation to Hq, thus 
- at least for low energies and far from the crossings - the 
eigenstates of the full spin Hamiltonian are close to that 
of Hq. As a consequence of the labelling convention in- 
troduced after Eq. ||SJ|, a given adiabatic eigenstate |i?„) 
before and after the avoided crossing of the levels En and 
En±i corresponds to two diflFerent states \m) ^ |m'), see 
Fig. n In other words, if the population corresponding 
to a certain adiabatic state does not change while pass- 
ing a crossing, the expectation value of Sz and thus the 
magnetization does change. 

The fact that the level splittings at the avoided cross- 
ings can differ by 12 orders of magnitude, raises an addi- 
tional difficulty when using Eq. ||7J) to calculate the dy- 
namics, because the derivatives also change in a similarly 
wide range. Therefore it turned out that a combination 
of Eqs. I0J and iQ leads to the most efficient method: 
when some populated levels are too close to each other, 
Eq. (UI) is no longer able to provide the required accuracy, 
therefore we change the basis and use Eq. Q for a short 
time interval after which it will be safe again to work with 
Eq. {T)). Thus the control parameter defining the stepsize 
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External magnetic field (T) 

FIG. 1: The level scheme of the Hamiltonian as a func- 
tion of the external magnetic field B. The energy levels 
corresponding to the exact eigenstates have "zigzag" 

form, while the approximate eigenenergies corresponding to 
the eigenstates of Sz are almost straight lines (with periodi- 
cally changing black and grey sections) for low energies. (A 
few examples are labeled in the figure.) The inset shows the 
minimal distance between levels as a function of B, the hori- 
zontal position of the crosses coincide the avoided level cross- 
ings shown in the main part of the figure. 



needed to have the required accuracy for such a long cal- 
culation is essentially the minimal distance between the 
populated levels. 

The dynamical equation jT)) indicates that the usual 
approach of treating the problem as a sequence of two- 
level transitions (each taking place at the corresponding 
level crossing) may provide an accurate approximation 
to the time evolution. In this framework the dynam- 
ics of the states corresponding to the anticrossing levels 
is governed by a 2 x 2 Hamiltonian resulting from the 
reductionii21SiSSiSL22iSa of the complete Hs to the rel- 
evant level pairs. Additionally, in this approach it is usu- 
ally assumed that the time dependence of the diagonal 
elements of the reduced Hamiltonian is linear, while the 
offdiagonal ones are constants: 



(9) 



With these assumptions each avoided level crossing is 
identical to the Landau-Zener-Stiickelberg^^'-'^^'^^ (LZS) 
model, which has an analytical solution yielding the tran- 
sition probability Plzs = 1 ~ exp(— 7rA^/257) in the 
long-time limit. Note that small values of Plzs means 
no appreciable change either in the population of the 
eigenstates of Sz, or the magnetization (but almost com- 
plete exchange of the populations of the adiabatic states); 
Plzs ~ 1 is observable as a step in the magnetization, 
while the populations of the adiabatic levels are practi- 
cally unchanged. Corrections to the LZS model originat- 
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FIG. 2: The expectation value of (Sz) (solid lines) as a func- 
tion of the external magnetic field B for different sweep rates. 
The exact results are compared with the predictions of the 
LZS approximation (dashed lines). The pulse shape corre- 
sponding to this figure is given by Eq. llObt . Above 5.7 T 
rapid oscillations appear, see the text for more details. 



ing from dipolar interactions have been investigated in 
ReL ^|. It is important to emphasize that Plzs de- 
pends on the ratio A/^/Q, i.e., on a single parameter 
(which, in appropriate dimensionless units, is simply the 
level splitting). In the next section we show that the dy- 
namics in the whole spin Hilbert-space can no longer be 
described by a single parameter, and consequently the 
exact transition probabilities can be significantly differ- 
ent from Plzs- 



III. UNITARY TIME EVOLUTION 

In this section we calculate the unitary dynamics de- 
scribed by Eq. |0J. Initially the external magnetic field 
is zero, then it raises to its maximal value of Bmaxj 
i.e., B{t) — /(t)-Bmax. Here we consider three analytical 
shapes of the fimction f{t): a linear, a sine and a tangent 
hyperbolical pulse, see the inset in Fig. |3| We construct 
these pulses in such a way, that the maximal external 
magnetic field rate w = max^dB / dt) = Bmax Taax{df /dt) 
is the same and falls in the kT/s range: 



hit) 


~ wt 


h{t) 


= sin 


hit) 


1 
2 



wt 



tanh( 



2wt - (5, 

R ' 



(10a) 
(10b) 

(10c) 



where the shift 5 in /3(t) has to be chosen such that at 
t = the external magnetic field is negligible. 

The initial state at the beginning of the calculation 
(t = 0) is the lowest energy eigenstate that later crosses 
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FIG. 3: The expectation value of (Sz) as a function of the 
time dependent external magnetic field B{t), which is shown 
in the inset, for three different pulse shapes with w = 2 kT/s. 



other adiabatic states, i.e., j^*) (0) — \Ei) w |m = —10). 
This means that we follow the lowest increasing curve 
on the level scheme shown in Fig. ^ and the energy lev- 
els that cross this line correspond to decreasing energies 
and thus do not meet any other levels later. Similarly, 
if initially the ground state \Eq) was populated, it would 
not give any contribution to the steps in the magnetiza- 
tion curve, it simply leads, to a very good approximation, 
to an additional constant. (Which is the reason for our 
choice of the initial state.) 

The expectation value (Sz) as a function of the ex- 
ternal magnetic field B is shown in Fig. [3 for the sine 
pulse Ijl0b|l and different values of the maximal sweep 
rate w. The steps seen in this figure are very similar to 
the experimental curves, but differ from the result that 
can be obtained by using the LZS theory (also plotted 
in Fig. O. Faster sweep rates mean smaller transition 
probabilities between the eigenstates of Sz ■ Although the 
exact dynamics is different from the LZS result, in the 
investigated sweep rate range we found that (Sz) scales 
with the sweep rate almost exactly the same way as one 
could deduce from Plzs- Additionally, Fig. [S] also shows 
that since the states |m) are not exact eigenstates of the 
complete spin Hamiltonian Hs, there are rapid oscilla- 
tions in (Sz) for higher external fields, which are clear 
indications of the Bohr oscillations corresponding to dif- 
ferent eigenenergies of Hs- We will find that if we take 
relaxation effects into account (Sec. IIV|I . these oscilla- 
tions disappear on a very short timescale. 

Comparison between the results for (Sz) for different 
functional shapes of the external magnetic field B(t) is 
shown in Fig.|31in case of a maximal sweep rate of w = 2 
kT/s. Note that the difference is not too large; in this 
sweep rate range it is not the functional time-dependence 
of B{t) that determines the heights of the steps seen in 
the magnetization curve, but rather its time derivative at 



the avoided level crossings. In other words, the approx- 
imation of a linearly increasing field B around a certain 
transition point is sufficient to accurately describe the 
dynamics at that transition. 

The population of the different eigenstates of Sz and 
the adiabatic states \En{t)) are shown in Fig. 21 as a func- 
tion of B for the representative example of a pulse with 
linear shape and w — 1 kT/s. As we can see, at the 
beginning of the time evolution, when the splitting of 
the adiabatic levels are very small and the magnetiza- 
tion is almost constant, the population of the |m = —10) 
state does not change significantly either. Around B ~ 5 
T the tunneling probability between different states |m) 
and |m') becomes appreciable, leading to the steps seen 
in Figs. [3 and 13 Note that the population of the adi- 
abatic levels show an opposite behavior: initially prac- 
tically all the population of the lower adiabatic level is 
transferred to the higher one at the avoided crossings. On 
the other hand, for larger external field values a nonzero 
population remains on the lower adiabatic level, leading 
to a noticeable change of the magnetization. The rea- 
son for the rapid oscillations in the populations of |m) 
seen in Figs. I2I4I is that at that high field strength Sz 
and Hs do not commute, therefore the states \m) are 
not exact eigenstates of the complete spin Hamiltonian. 
These oscillations disappear when we include relaxation, 
see Sec.Hvl 

If we restrict ourselves to the LZS model, then starting 
from the ground state it is sufficient to find the first value 
of B at which this adiabatic level anticrosses the next 
one, calculate the LZS parameter A/^/Ti and use Plzs 
to obtain the population of the two relevant adiabatic 
levels after the transition, and repeat this process until 
the end of the time evolution. This approach includes 
the slow change of the levels \E„{t)) , causing a slight, 
continuous increase of the magnetization as a function of 
B, but this effect is difficult to see in Figs. El and 13 the 
steps dominate the behavior of (Sz)- However, the LZS 
result obtained in this way is quantitatively different from 
the exact {Sz)it) curve that was calculated by taking all 
the 21 levels into account (see Fig. O- The position of 
the steps (determined by the avoided crossings) are the 
same, but their heights are different, and this difference 
can be as large as 30%. 

In order to understand the physical reason for this ef- 
fect, we have to investigate the relation between Plzs at 
a certain anticrossing and the relevant transition proba- 
bility resulting from the present calculation. The first im- 
portant point to take into account is that, to a very good 
approximation, the transitions seen in Fig. ^ take place 
between two neighboring adiabatic levels. For sweep 
rates in the kT/s range the characteristic time of the 
transitions'^^ at the avoided level crossings neither over- 
lap nor influence each other. Let us now concentrate 
on a single step, as an example we focus on the vicinity 
of i? = 5.7 T (see Fig. jSJ, where the anticrossing levels 
Ell and £'12 correspond to the states \Eii) « |m = —10) 
(|m = -l)) and \Ei2) « |m = -1) (|m = -10)) before 
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a) 




FIG. 4: /Reduced resolution for this preprint version./ The 
population of the levels corresponding to the states \En) [a)] 
and |m) [b)] as a function of the external magnetic field for a 
linear pulse with w = 1 kT/s. For the sake of simplicity in b) 
we show only populations larger than 0.05. 

(after) the transition. 

The most remarkable point concerning Fig. is that it 
was obtained by assuming that only the adiabatic levels 
and \Ei2) play a role in the transition, the dynamics 
has been reduced to these two relevant levels. Calculat- 
ing the same transition by taking all the 21 levels into 
account shows that this approximation estimates the ex- 
act dynamics with high accuracy. More generally, it is 
always possible to perform a similar reduction around a 
given avoided level crossing. Thus we obtain a method 
where a sequence of effective two-level transitions can 
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FIG. 5: The population of the level En around the transition 
— > |i5i2) for different sweep rates. We compare the 
exact results with those from the LZS model. The external 
field pulse corresponding to this figure is given by Eq. (llUal . 



describe the time evolution. The results obtained in this 
way are practically the same as those of the exact calcu- 
lation, and considering the numerical costs, this is a very 
effective method. 

However, the time dependence of the expectation value 
of (Sz) obtained in this way is still different from the LZS 
result, despite the fact that the LZS theory is also based 
on a two-level approximation. Note that all the curves 
shown in Fig. were calculated using the same numer- 
ical method, the pulse shape (/i), the initial state and 
the level splitting were also the same for all the calculated 
sweep rates: the only difference was the time dependence 
of the energy levels and their coupling. (In fact, as it is 
clear from Eqs. Q and JSJ, it is only the energy 
difference Ei2{t) — Eii{t) that plays a role here.) These 
are multilevel effects: the time dependence of the 2x2 
Hamiltonian obtained by the reduction of Hs to the rele- 
vant level pair is affected by all the other levels, similarly 
to a renormalization effect. The influence of the states 
not taking part in the transition results in a time de- 
pendence of the parameters of the reduced Hamiltonian 
which is slightly different from the LZS model described 
by Eq. J^. That is, the single parameter A/-\/?2 is not 
enough to describe these transitions, or in other words, 
we have time dependent factors in the LZS matrix ele- 
ments A(i),f2(t). That is emphasized by calculating the 
long time limit transition probability for the same tran- 
sition shown in Fig. [S] with different parameters in the 
Hamiltonian (PJ in such a way, that the sweep rate and 
the level splitting are fixed. As an example, we consider 
the level splitting A at this transition as a function of 
only two parameters, namely C and K] we do not change 
D, F and E in Hs- First we calculate A(C, K) using the 
parameters given in Sec.lTTland obtain Aq = 1.47 x 10~^ 
(in Kelvin). Then we solve the equation A(C, if) = Aq 
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5.0x10' 1.0x10' 1.5x10' 2.0x10' 

FIG. 6: The final population of the level E\i after the tran- 
sition l-Eii) ^ \Ei2) around B = 5.7 T as a function of the 
parameter C in Eq. ||3J. The level splitting min(i5i2 — 
is kept fixed, leading to the constant final population in the 
LZS model indicated by the solid horizontal line. The mag- 
netic field pulse used is given by Eq. llOal with w = 1 kT/s. 

to find the parameter pairs {C, K} that correspond to the 
same level splitting as it was initially. It turns out that 
there are several disconnected lines in the C — K plane 
along which A{C,K) = Aq, and consequently Plzs is 
constant for a given sweep rate. But as we can see in 
Fig. IHI the calculated non-LZS transition probabilities, 
long after the transition took place, strongly depend on 
the parameter C in Eq. ||2Jl: they range from 5% up to 
65%, and can be both higher and lower than Plzs- Dif- 
ferent symbols in Fig.|Slcorrespond to different lines along 
which A{C,K) — Aq. In any case, we can conclude that 
the final transition probability at a given avoided cross- 
ing is also influenced by the levels that do not take part 
in the relevant transition, and can therefore not be de- 
scribed within the framework of the LZS model. 



IV. RELAXATION EFFECTS 

So far we considered unitary time evolution, i.e., the 
Hamiltonian ^ governed the dynamics. However, there 
are interactions with e.g. phonons that are not included 
in Hs- These additional degrees of freedom can be con- 
sidered as the environment of the spin system, and 
any realistic description should take their influence into 
account Additionally, as we shall see in this sec- 
tion, the rapid oscillations seen in Figs. 12141 disappear on 
a very short time scale when relaxation influences the 
dynamics. 

Since phase relaxation is usually much faster than en- 
ergy exchange between the investigated quantum system 
and its environment, we concentrate on this kind of deco- 
herence, and assume a Lindblad-type-i^ dynamical equa- 
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FIG. 7: The effect of phase relaxation on the dynamics of the 
magnetization ((Sz)) around the transition |-Bii) \Ei2). 
We assumed a linear pulse given by Eq. HlOajl with maximal 
sweep rate w = 1 kT/s in part a) and 2 kT/s in part b). 
Note that the origin of the time axis has been shifted, t — 
corresponds to the time when E12 — En is minimal. 

tion: 

= [Hs, + I {2S,gS. - S'^g - gS',) . (11) 

The second term in this master equation will not change 
either the magnetization of the sample nor the expecta- 
tion value of Hq. It leads to the gradual disappearance 
of the non-diagonal elements of g in the eigenbasis of 5*2 
without changing the populations . The result of this 
kind of relaxation is quite different during the transitions 
at the avoided level crossings and between them. This 
difference is clearly seen if we consider sweep rates of the 
order of kT/s, when the characteristic transition times 
are 10-6-10-'' s, while the spin system spends about 
IQ-'^ s between two level crossings. These time scales 
should be compared to 7-^ , which is in the range of IQ-^- 
IQ-^ s for low temperatures (T w 2K) at which several 
experiments were performed. Thus, between two cross- 
ings, relaxation has enough time to destroy the phase 
information almost completely. 

In other words, if after an avoided crossing the system 
is in a pure quantum mechanical state \(f>), the initial 
condition at the next crossing can be considered as the 
right hand side of the following scheme: 

I'/') {(t'l = X! 1""^ ^"^1 X! 1"^^ ^"^1 ^'^'^^ 

nm m 

Note that the Hamiltonian part of the time evolution 
slightly modifies the process above, but the final con- 
sequence, i.e., that dephasing changes the initial condi- 
tions for the transitions at the avoided level crossings is 
still valid. Decoherence as described by Eq. H12|) plays 
an important role in the physical mechanism responsible 
for the observed form of hysteresis loops (see Ref. (S^l 
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for experimental results on systems of low-spin magnetic 
molecules). We found that even for such fast external 
magnetic field sweep rates as a few kT/s, quantum me- 
chanical interference does not play an important role in 
transitions at consecutive anticrossings, by the time the 
system reaches the next transition point, the phase infor- 
mation has already flown into the environment. For sim- 
ilar reasons the formation of closed hysteresis curves can- 
not be influenced either by quantum interference. How- 
ever, even faster sweep rates, or faster return to a certain 
transition by superimposing an oscillating magnetic field 
on a constant one around a certain crossing, may give 
rise to quantum mechanical interference effects. We pre- 
dict that such interference phenomena will show up for 
magnetic field sweep rates in the MT/s range. They will 
be seen as a very strong dependence of the form of the 
hysteresis loops on the sweep rate due to the extreme sen- 
sitivity of the process to the relative phases of the states 
that take part in a transition. 

The phenomena expected around transition regions are 
more interesting, as they are consequences of the inter- 
play between coherent effects and dephasing: the charac- 
teristic time of the transitions are comparable, but usu- 
ally shorter than the dephasing time defined by 7""'^. The 
general consequence of the second term in the master 
equation is that it decreases the transition probabil- 
ity in the {Iti)} basis, i.e., when relaxation is present, it 
makes the steps in the magnetization smaller (see Fig.[71l. 
The larger the value of 7, the stronger this effects is, 
which underlines again the importance of the relative 
phase in this physical system. In the previous section 
it was shown that phases gained by the adiabatic states 
at an anticrossing can remarkably modify the transition 
probability, now we see that loss in phase information has 
also strong effects. Note that in contrast with the coher- 
ent case (see Sec. the final transition probabilities 
in the presence of dephasing do not scale with the sweep 
rate similarly to the LZS case. For larger sweep rates 
the system spends less time in the transition region, and 
consequently decoherence less strongly modifies the final 
transition probability. That is, the larger the magnetic 
field sweep rate, the more similar the dynamics around 
a transition becomes to the case without dephasing. Ad- 
ditionally, we want to point out that the coherent os- 
cillations, the consequences of which has been seen in 
Figs. I2I4I are strongly damped even for weak dephasing. 

However, the statement that the final transition prob- 



abilities are modified by the presence of all the levels and 
thus cannot be accurately described within the frame- 
work of a LZS model (even if we include relaxation 
as well) is also true in the present case. Additionally, 
let us emphasize that the parameter region discussed 
here is d ifferen t from the strongly damped one studied 
m Refs. '28*35|, where incoherent tunneling can give a 
proper description. For external field sweep rates in the 
kT/s range, neither coherent time evolution nor relax- 
ation dominates, which leads to an interesting interplay 
between these two qualitatively different processes. 
V. CONCLUSIONS 

We studied the time evolution of the spin degrees 
of freedom in molecular nanomagnets, with a focus on 
the molecule Mni2-Ac in the presence of time-dependent 
magnetic field. Using an appropriate 'exact numerical' 
method, we followed the time evolution from zero exter- 
nal magnetic field until saturation of the magnetization is 
reached. We found that for sweep rates in the kT/s range, 
steps in the magnetization originate from two-level tran- 
sitions which cannot be described within the framework 
of the Landau-Zener-Stiickelberg (LZS) model. This ob- 
servation led us to the introduction an efficient and ac- 
curate approximation based on two-level non-LZS transi- 
tions. This method introduces the possibility of perform- 
ing long term dynamical calculations that can directly 
be related to experiments. We also demonstrated that 
the sweep rate range of kT/s is special in the sense that 
for realistic relaxation times, there are observable conse- 
quences of the competition between coherent dynamics 
and decoherence that modify the heights and widths of 
the magnetization steps. 



Acknowledgments 

We thank V. V. Moshchalkov and J. Vanacken for 
the fruitful discussions and O. Kalman for her valuable 
comments. This work was supported by the Flemish- 
Hungarian Bilateral Programme, the Brazilian Council 
for Research (CNPq), the Flemish Science Foundation 
(FWO-Vl), the Belgian Science Policy and the Hungar- 
ian Scientific Research Fund (OTKA) under Contracts 
Nos. T48888, D46043, M36803, M045596. 



^ D. Gatteschi, R. Sessoli, and J. Villain, Molecular Nano- 
magnets (Oxford University Press, 2006). 

^ J. R. Friedman, M. P. Sarachik, J. Tejada, and R. Ziolo, 
Phys. Rev. Lett. 76, 3830 (1996). 

^ K. M. Mertes, Y. Suzuki, M. P. Sarachik, Y. Paltiel, 
H. Shtrikman, E. Zeldov, E. Rumberger, D. N. Hendrick- 
son, and G. Christou, Phys. Rev. Lett. 87, 227205 (2001). 

* W. Wernsdorfer, M. Murugesu, and G. Christou, Phys. 



Rev. Lett. 96, 057208 (2006). 

^ J. Tejada, E. M. Chudnovsky, J. M. Hernandez, and 
R. Amigo, Appl. Phys. Lett. 84, 2373 (2004). 

^ J. Vanacken, S. Stroobants, M. Malfait, V. V. 
Moshchalkov, M. Jordi, J. Tejada, R. Amigo, E. M. Chud- 
novsky, and D. A. Garanin, Phys. Rev. B 70, 220401 
(2004). 

^ A. Hernandez-Minguez, M. Jordi, R. Amigo, A. Garcia- 



8 



Santiago, J. M. Hernandez, and J. Tejada, Europhys. Lett. 

69, 270 (2005). 
* R. M. Dickc, Phys. Rev. 93, 439 (1954). 
" M. G. Benedict, A. M. Ermolaev, V. A. Malyshev, I. V. 

Sokolov, and E. D. Trifonov, Superradiance (lOP, Bristol, 

1996). 

^° E. M. Chudnovsky and D. A. Garanin, Phys. Rev. Lett. 

89, 157201 (2002). 
" V. K. Henner and I. V. Kaganov, Phys. Rev. B 68, 144420 

(2003). 

C. L. Joseph, C. Calero, and E. M. Chudnovsky, Phys. 
Rev. B 70, 174416 (2004). 

M. G. Benedict, P. Foldi, and F. M. Peeters, Phys. Rev. B 
72, 214430 (2005). 
" Y. Suzuki, M. P. Sarachik, E. M. Chudnovsky S. McHugh, 
R. Conzalez-Rubio, N. Avraham, Y. Myasoedov, E. Zel- 
dov, H. Shtrikman, N. E. Chakov, et al., Phys. Rev. Lett. 
95, 147201 (2005). 

A. Hernandez-Mmguez, J. M. Hernandez, F. Macia, 
A. Carcia-Santiago, J. Tejada, and P. V. Santos, Phys. 
Rev. Lett. 95, 217205 (2005). 
^'^ M. N. Leuenberger and D. Loss, Nature (London) 410, 789 
(2001). 

L. D. Landau, Phys. Z. Sowjetunion 2, 46 (1932). 
^® C. Zener, Proc. Roy. Soc. London, Ser. A 137, 696 (1932). 
1^ E. C. G. Stiickelberg, Helv. Phys. Acta 5, 369 (1932). 
^° I. Mirebeau, M. Hennion, H. Casalta, H. Andres, H. U. 



Ciidcl, A. V. Irodova, and A. Caneschi, Phys. Rev. Lett. 

83, 628 (1999). 

A. L. Barra, D. Gatteschi, and R. Sessoh, Phys. Rev. B 
56, 8192 (1997). 

S. Hill, J. A. A. J. Perenboom, N. S. Dalai, T. Hathaway, 
T. Stalcup, and J. S. Brooks, Phys. Rev. Lett. 80, 2453 

(1998) . 

S. Hill, R. S. Edwards, S. I. Jones, N. S. Dalai, and J. M. 

North, Phys. Rev. Lett. 90, 217204 (2003). 
'^^ J. H. van Vleck, Phys. Rev. 33, 467 (1929). 

J. des Cloizeaux, Nucl. Phys. 20, 321 (1960). 
2® D. J. Klein, J. Chem. Phys. 61, 786 (1974). 

D. A. Garanin, J. Phys. A 24, L61 (1991). 
^® M. N. Leuenberger and D. Loss, Phys. Rev. B. 61, 1286 

(2000). 

S.-K. Yoo and C.-S. Park, Phys. Rev. B 71, 012409 (2005). 
^° J. Liu, B. Wu, L. Fu, R. B. Diener, and Q. Niu, Phys. Rev. 

B 65, 224401 (2002). 
3^ N. V. Vitanov, Phys. Rev. A 59, 988 (1999). 

I. Rousochatzakis and M. Luban, Phys. Rev. B 72, 134424 

(2005). 

G. Lindblad, Commun. math. Phys. 48, 119 (1976). 

M. Vogelsberger and D. A. Garanin, Phys. Rev. B 73, 

092412 (2006). 

M. N. Leuenberger and D. Loss, Europhys. Lett. 45, 692 

(1999) . 



